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■ In this paper we study a new family of sinc-like functions, defined on an interval of finite width. 

' These functions, which we call "little sine", are orthogonal and share many of the properties of the 

sine functions. We show that the little sine functions supplemented with a variational approach 
enable one to obtain accurate results for a variety of problems. We apply them to the interpolation 
of functions on finite domain and to the solution of the Schrodinger equation, and compare the 
performance of present approach with others. 
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I. INTRODUCTION 



The main goal of this paper is to introduce a new set of orthogonal functions, hereafter called "little sine functions" 
(LSF), and show that they can be used to solve a wide class of problems, such as function interpolation, the numerical 
solution of differential equations, including the Schrodinger equation. Actually the LSF, which arc defined on a finite 
domain, were first derived by D. Baye [lj in the framework of generalized meshes, but the relation with the usual 
sine functions was not recognized. As a matter of fact, in recent years Baye and collaborators have used these and 
£7^" other functions for the solution of the Schrodinger equation with several different potentials, and they have produced 
' accurate numerical results both for the energies and wave functions [HIE El- 
• i-h , On the other hand, sine collocation methods have been applied to a larger class of problems, which include those 
' mentioned above (see for example Q), and they have been put on firm mathematical grounds Q. Examples of 
applications of the sine functions can be found in references 0, 0, H, EJ El [pi , among others. Recently, one of the 
authors has shown that the sine collocation methods can be optimized by using a variational approach based on 
the Principle of Minimal Sensitivity |l3j . This optimization allows one to obtain the highest precision with a given 
number of grid points, and can be particularly valuable in problems that require intense numerical calculation. 

We feel that it is important to establish the relation of the LSF with the sine functions at least for two reasons: 
first, because it will be possible to generalize the variational approach 01 to the LSF, thus improving the convergence 
behavior of the method, and second because it offers an interesting link to generalized mesh methods and to open 
new areas of application of the latter. 

This paper is organized as follows: in Section [n] we describe the general properties of the usual sine functions, 
defined on the real line. In Section [illl we derive an expression for the LSF and discuss their properties. In Section llVl 
we solve the Schrodinger equation with two different potentials and compare numerical results obtained with the usual 
sine functions and LSF. Finally, in0we draw our conclusions. 

II. SINC FUNCTIONS 

In what follows we outline some of the basic properties of the generalized sine functions, defined as 

sm(n/h(x-kh)) 
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where k — 0, ±1, ±2, .... The sine function for a given value of the index k is peaked at x k = kh, where it equals 
unity, and vanish at the other points Xj = jh, with j ^ k and j = 0, ±1, ±2, .... 
From the integral representation 



S k (h,x) = A / e ±*(*-*fc)* rff , (2) 
we easily derive the normalization factor 

/■+O0 

Ii = Sk{h,x) dx = h , (3) 



and the orthogonality property 

" + OG 

Sk(h,x) Si(h,x) dx = h Ski- (4) 

> 

It is worth noticing that cq. J5J defines a Dirac delta function in the limit h — > 0. 

A function/ (x) analytic on a rectangular strip centered on the real axis can be approximated in terms of sine 
functions as 



/(*)« E -f( fc/i ) > ( 5 ) 



4-oc 



k— — oo 

which together with the normalization factor can be used to approximate the definite integral 

r*+oo 



/-|-oo _ _ 
f(x) dx ^ h f( k > h )- (6) 

- OO / 



fc— — oo 

It is not difficult to derive simple expressions of the derivatives of sine functions in terms of the same sine functions: 



—S k (h,x) = ^ c\l> S k (h,x) , (7) 



dx 

1 — — OC 



where 

„(i) - 



+oo i if k = I 

^ S l (h,x)-S k (h,x)dx=l i f/ M/ (8) 



For the second derivative we have 



j2 00 

^ s S h (h,x)= E 4? , (9) 

/— — oo 



where 



d^ . if k — I 



4 2) = / ^s k (h,x) dx = _ 2 ( 3 _\>-r ., , (io) 



'W (k-l) 

General expressions for higher order derivatives are also available 

J?r) — (-l)'~ fc W iNfe+l 2r ' „2fc /. ,\2fc 



W7 



(if 



(2r+l) = (-1) 



II 

with r = 1,2,.... 
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FIG. 1: LSF for N = 20 and different values of k (color online). 



III. LITTLE SINC FUNCTIONS 



We will now derive the new set of LSF. We consider the orthonormal basis of the wave functions of a particle in a 
box with infinite walls located at x = ±L: 

V>„(z) = sin(^(a; + L)) (15) 



and define 



n ( zin( {2N+1)7r(x - y) ) coc ( (2N+1)n{x+y) ) 

M*, y)-c N Y: i> n (x)Mv) - if V ( ( 4L A J (-i) w V ( " A J 



(16) 



4L I 



where Cat is a constant. 

Because of the completeness of the basis {ip n (x)} we have 

6 N (x,y) . ■ , 
J im — = °i x ~ V) ■ ( 17 ) 

For reasons that will soon become clear, we set Cjv = ^ and select even values of N. To simplify the notation 
h = 2L/N will denote the grid spacing, and y k = ^ = kh, with = -N/2 + 1, -iV/2 + 2, . . . , N/2 - 1, the grid 
points. 

We then define a set of N - 1 LSF 

s (h N x) = 1 / sin((l + ^) j{x-kK)) cos((l + ^) l(x + kh)) } 

27V \ sin^^-fcft)) cos^x + kh)) j' 

Fig- ID shows 5 of these functions for TV = 20. 

It is not difficult to prove that the are orthogonal 

L 

Sk(h,N,x)sj(h,N,x)dx = hdkj ■ (19) 

J-L 

and satisfy 

s k (h,N, yj )=S k j , (20) 
properties that are also shared by the sine functions. 
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FIG. 2: Comparison between the sine functions and little sine functions corresponding to N — 40 (color online). 



Therefore, it is not surprising that the LSF become the standard sine functions when N —> oo and h is held constant 
in eq. (JTSJ: 



/, »t \ &ixi(ir(x — kh)/h) „ ., . 
hm s k (h,N,x) = vv lh J!' = S k (h,x) 
N^oo TT{x—kh)/h) 



(21) 



This property justifies the name of little sine functions. 

Fig. [3 shows the LSF for N = 40 and L = 1 together with the corresponding sine functions. Differences between 
both kind of functions are appreciable only in the right plot, corresponding to k = 19. In this case the LSF is slightly 
larger than unity at the peak and its oscillations die out faster. 

The LSF share some properties with the sine functions; for example, we can approximate a function f(x) on the 
interval (-L, L) as 



N/2-1 

/(*)« E f{xk)a k (h,N,x), 

k=-N/2+l 

where Xk = kh. Similarly one can express the derivatives of LSF in terms of LSF as 
ds k (h,N,x) r-\ ds k (h,N,x) 



dx 

d 2 Sk(h, N, x) 
dx 2 



E 



d.r 



E 



d 2 s k (h,N,x) 



dx 



Sj (h,N,x)=J2 c >3 s 3 {h,N,x) 

3 3 



(22) 

(23) 
(24) 



(r) 

where the cL are the counterpart of the coefficients shown above for the sine functions. 



An explicit calculation yields 



c (1) 
c kj 



and 



4i 



tan 



cot 

V ' 4L 1 



U - k)TT 

2N 



tan 



(j + k)ir 
2N 



,(2) _ 



JJ 



,(2) _ 



24L 2 



cos (f) cos (^) 



8L2 cos 2(_^ (j+fc) ) sin 2 ( _^ C? _ fc)) 



(25) 
(26) 

(27) 
(28) 



Because of eq. (12111 these matrices reduce to the usual sine expressions given in the preceding section when N — > oo. 
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It is straightforward to generalize from cq. I|18f) . defined in (— L, L), to an arbitrary interval (a, b): 



s k (h,N,x) = s k h,N, 



2L 



a + b 



In this case the points of the grid are given by 



x k = (b-a) 



N 



(29) 



(30) 



In order to apply sine collocation on a finite domain one commonly maps the real line onto a finite interval || by 
means of the conformal transformation 



(z) = log 



This map carries the eye-shaped region 

De = <j z 

into the infinite strip 



arg 



77 

<d<- 



Ds = = u + iv : \v\ < d < — j 



(31) 



(32) 



(33) 



Under the inverse transformation, z = (f>^ 1 (w) the points of the uniform grid on the real axis, given by u k = kh, 
are mapped onto the non uniform grid defined by the points x k = (a + be kh )/(l 
In this case the sine functions are mapped onto 



Sk(h, x) 



sin (ir((fi(x) — kh)/h) 



n((j)(x) - kh)/h ' 

which equals unity at x = x k . Consequently, it is possible to approximate a function in the interval (a, b) as 

+N 

f(x)*7(x)= m)s k (h, X ) . 

We can test the performance of the LSF on an example selected from Ref. : 

f(x) = 2x 2 +x- 3x 3 



(34) 



(35) 



(36) 



where < x < 1. Fig. [3] compares the logarithmic error A(x) = log 10 \f(x) — /(x)| for both kind of functions. The 
solid curve corresponds to 21 LSF, whereas the other three curves correspond to the same number of conformally 
mapped sine functions with spacing h = 1/4, 1/2, and 1 respectively. The LSF produce smaller errors and offer the 
advantage of a uniform grid. 

Stenger originally introduced sine methods for the numerical solution of differential equations 0. Sinc-Galerkin 
and sine collocation methods arc particularly useful in dealing with these problems since they converge exponentially 
even in the presence of boundary singularities. It is not our purpose to generalize all the known mathematical results 
from the sine functions to the LSF; however, we assume that both kind of functions share similar properties and 
simply compare our LSF results with those provided by the conformally mapped sine functions. 

We consider the example 4.1 of Ref. @; that is to say, the inhomogeneous differential equation 1 



u"(x) + u'(x) + u(x) 



(I) 2 (■ 



2x 3 - 29x 2 + 62x + 38) 



(37) 



1 Notice a typo in 0, where 69a; must read 62a;. 
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with the boundary conditions u(—l) = it(4) = 0. The exact solution to this equation is 



Uexact(x) = ( ^ ) {x + if {x - 4) 2 



We look for a numerical solution in terms of our LSF as 

JV/2-1 



ii(i)«mi s (i)= ^2 Uk 8k{h, N,x - 3/2) 

k=-N/2+l 



Fig. 0| shows global and local errors defined respectively as 

r 4 

E G (N) = log 10 



log 



10 



{Uexact{x) ~ Ui s (x)) 2 dx 
J —1 I — Ar/o 



fe=-AT/2+l 

S L (Af) = log 10 |(u e:C act(3/2)-Ki s (3/2)| . 



(38) 



(39) 



(40) 
(41) 



Notice that at large N the errors appear to decay exponentially. Unlike Ref. Q no domain decomposition was 
required to solve this problem. 



IV. THE SCHRODINGER EQUATION 

As mentioned in the Introduction, sine collocation methods have been used to obtain accurate numerical solutions 
to the Schrodinger equation. In this section we wish to extend the variational approach and the results of Ref. 01 
to our LSF and generalize the method enclosing potentials with both bound and unbound states. To simplify this 
presentation we begin with potentials containing only bound states. 

A. Potentials containing only bound states 

We consider the one dimensional Schrodinger equation 

h 2 d 2 ^ n {x) 



+ V{x)i> n {x)=E n i> n {x) (42) 



5 I i I i I i I i I i I 

20 40 60 80 100 

N 



FIG. 4: Global and local errors I4U1 and I4U . respectively, for the solution of eq. 1371 in terms of N — 1 LSF. 




FIG. 5: Comparison between the exact solution utrueix) (solid line) and the approximations obtained using LSF with N = 
4, 6, 10 (color online). 



on an interval (a, 6), which can be either finite or infinite. The wave functions ip n {x) obey the boundary conditions 
"0n( a ) = ^nip) = which grant that there will be only bound states. 

If we express the wave functions in terms of our LSF. then equations 1221) . (|24|l and (|28|l allow one to derive the 
following matrix representation of the hamiltonian operator 

H k i = -^-c^> +5 kl V(kh) . (43) 
2m 

This equation is similar to eq.(8) of Ref. except for the form of the matrix c^ 2 ). Once the spacing h of the 
grid is fixed, then the N x N matrix Hjy for a manifold of N LSF can be diagonalizcd. In this way one obtains 
an approximation to the lower part of the spectrum, consisting of the first N eigenvalues and wave functions of the 
Schrodinger equation (|42|l . 

Going back to eq. I|43|) we observe that the precision of the approximate results obtained by diagonalization of the 
hamiltonian matrix depends crucially on both the number of LSF as well as on the grid spacing. In fact, altough 
a small spacing can help to increase the resolution, if the number of the sine functions is not large enough, the 
approximation will not be able to grasp the natural scale of the problem and the overall precision will be poor. On 
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the other hand, a large spacing will certainly provide poor results, because the details of the problem will not be 
resolved. On these grounds it is easy to convince oneself that it is likely to exist an optimal spacing, for a given 
number of functions. Finding this optimal spacing will allow one to reach sufficiently accurate results with a relatively 
small number of grid points. 

It was found earlier [l2^ that an optimal spacing can be obtained by straightforward application of the principle 
of minimal sensitivity (PMS) to the trace of the N x N hamiltonian matrix 7jv(/i) = Tr [Hjv]- In fact, given that 
the trace of a hamiltonian is invariant under unitary transformations, and that in the limit h — > it will therefore be 
independent of h, then the optimal spacing may be given by the solution of the equation 2 

±T N {h) = Q. (44) 

Therefore the optimal value of h is found by numerically solving a single algebraic equation, a modest computational 
task; thus, the interval length L appearing in LSF is treated as a variational parameter. 

To test the performance of our method we apply it to the first example in Ref. Q, i.e the harmonic oscillator with 
fi = 1, m =1/2 and u> = 2. Using a cartesian mesh, Baye and Heenen reported errors smaller than 10~ 3 for the 
first three eigenvalues and N = 10. On the other hand, our LSF approach with N = 10 (which corresponds to 9 
sine functions), and a grid spacing optimized according to eq. (|44|) . yields errors of —4.86 x 10 -6 , 1.2 x 10 -4 and 
— 1.6 x 10~ 3 for the same cases. 

The authors of Ref. Q also observe that for TV = 50, the high eigenvalues become very sensitive to the value of h, 
and that the variation of the 30th eigenvalue with respect to h presents a marked minimum around h = 0.35. It is 
remarkable that the PMS condition for N = 50 yields Lpms = 8.93, corresponding to hpMS = %Lpms/N as 0.357, 
which is extremely close to the value quoted by Baye and Heenen. It is worth noticing that while the optimal value 
of h quoted by Baye and Heenen is the result of an empirical observation, the almost identical value of h given by the 
PMS is just the numerical solution of the algebraic equation Q44(L which requires negligible computer time. 

As a second example of application of the PMS to the LSF collocation method we consider the anharmonic potential 

V(x) =x 2 +x 4 (45) 

and we assume that h = 1 and m = 1/2 in the Schrodinger equation. This example was also studied by Baye and 
Heenen Q who obtained the optimal values h = 0.55 for N = 10, and h = 0.2 for N = 50 using a cartesian mesh, 
which is somehow related to the LSF. Fig. [5] shows the error ^^j^ xac ^ — | as ct function of the parameter L 

for three different values of N (remember that the number of LSF in the expansion is N — 1). The plus symbols in 
the plot correspond to the predictions of the PMS condition, which generally fall close to the minimum of the curve, 
while the square symbols correspond to the solutions of a sort of empirical PMS condition, obtained by minimizing 
the modified trace 

N — 2 

f N {h)=T N {h)-^ r Y,V{x k ) . (46) 

k 

with respect to L. The better behavior of the modified PMS condition for just one state (as in Fig. |SJ) should not 
confuse the reader. It must be kept in mind that the PMS minimizes a sort of global error for all the states in the 
chosen manifold. In order to appreciate this point clearly, Fig. shows the global error a = J2n=o ~ En Xact ^\ 

as a function of L for the potential x 4 /4 and N = 20. Notice that in this case the PMS condition yields the minimal 
error. The exact energies j^ xact ^ were simply chosen to be those given by the method at higher order: En°\ 

The accuracy of present calculations is greater than that obtained earlier for the same problem Q, where the 
authors report errors of the order of 10 -5 and 10 -12 for N = 10 and N = 50, respectively. Fig. shows that the 
curves with different values of N overlap in the region of small L, which suggests that the approach may not be taking 
into account the large-L region correctly. 



B. Potentials containing bound and unbound states 

In what follows we try and show that the LSF are suitable for the treatment of potentials with bound and un- 
bound states. The application of the PMS condition to potentials with both discrete and continuous spectra is not 



2 Notice that such property was invoked earlier in Ref. 1 1 -41 when using a basis of harmonic oscillator wave functions depending on an 
arbitrary scale parameter 
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FIG. 6: Error in the ground-state energy of the potential 1451 1 using N = 10, 20 and 30, as a function of the parameter L. The 
plus and square symbols correspond to the PMS and modified PMS conditions, respectively (color online). 



O 10 3 r 




FIG. 7: Global error a = J^lLo I -En - E { „'\ for the potential V(x) = i 4 /4 (color online). 



straightforward. The potentials treated in the preceding section increase with the coordinate and therefore the matrix 
representation of those potentials increase with L. Since the matrix representation of the kinetic energy decreases 
with L then there is a minimum in the trace of the Hamiltonian matrix and the PMS gives a balance between the 
traces of the kinetic and potential energies. That minimum provides the natural length scale for the application of 
the method. If the potential energy tends to a finite constant value as the coordinate increases, then there may not 
be a minimum and the PMS will not yield the length scale for the application of present approach. 

In order to overcome that problem we substitute a potential V(x) into the Schrodingcr equation that behaves 
exactly as the original potential V(x) in the relevant coordinate region and increases to infinity at large distances. 
The error introduced by such potential substitute will be negligible if the difference between the original and substitute 
potentials is appreciable only where the wavefunctions are expected to be vanishing small. This substitution removes 
the continuous spectrum but should not affect the discrete one too much. The advantage is that we are thus able to 
apply the PMS condition exactly as in the preceding subsection. The potential substitute is introduced with the sole 
purpose of obtaining the length scale and we diagonalize the correct hamiltonian matrix. 

We will illustrate our procedure on the Morse potential already treated earlier by means of the Lagrange mesh 
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TABLE I: Errors 


rpapprox 

£n — -CWi 


- E e n xact for 


some states of the Morse potential. 


Powers of ten are indicated between 


square 


brackets. 












1 


n 


N 


Ref.[3] 


PMS 


hpMS 








20 


3.8[-7] 


-3.3[-10] 


0.223 






40 


< 


1.8[-20] 


0.134 




5 


20 


4.0[-3] 


2.1[-6] 


0.223 






40 


l.l[-9] 


2.5[-14] 


0.134 


1 





20 


4.1 [-7] 


-3.6[-10] 


0.223 






40 


< 1[~14] 


-1.6[-20] 


0.134 




5 


20 


4.4[-3] 


1.5[-6] 


0.223 






40 


l.l[-9] 


2.9[-14] 


0.134 


2 





20 


4.5[-7] 


-4.1[-10] 


0.223 






40 


< 1[-14] 


1.2[-20] 


0.134 




5 


20 


4.7[-3] 


2.0[-7] 


0.223 






40 


1.1 [-9] 


3.9[-14] 


0.134 



method 

V(r) = D le-Vr-r.) - 2e- a(r - r ^l , (47) 

where D = 0.10262, r e = 2. a = 0.72, 2m = 1836 and h = 1. In the case of states with nonzero angular momentum 
we should add the centrifugal potential > £ = 0,1, 2,.... 

The potential substitute F(r) is arbitrary; we can, for example, choose it to be the Taylor expansion of V(r) around 
a given point, and truncate it at a sufficiently large order to be accurate enough at small r, and at the same time to 
satisfy linv^oo V(r) = +oo. For example, in present particular case we can choose: 



£ — < 17! 



n\ dr r 

n=0 



(r - 10)" . (48) 



r=10 



Another difficulty to take into consideration is that our LSF arc defined in (-L, L) whereas the radial coordinate is 
defined in (0, +00). The obvious solution to this apparent problem would be using the form of Ref. where the left 
boundary of the interval is 0; on the other hand, such choice would be equivalent to a shift by a proper amount of the 
potential, thus bringing the boundary condition on the left point to coincide with the left point of the LSF. Despite 
its simplicity, this procedure is not optimal and generally does not provide the best results. We have found out that 
a more convenient strategy consists of keeping our LSF unchanged and shifting the coordinate by given amount as 
V(r) — > V(r + r), where r is typically close to the minimum of the potential (the same shift should be applied to the 
centrifugal energy when necessary). The PMS applied to the shifted hamiltonian will thus provide the optimal scale 
for the application of the LSF collocation method. 

One advantage of that procedure is that it maximizes the sampling of the classically allowed region, where the bound 
state wave functions exhibit marked nonzero contributions. Of course, in order to take into account the boundary 
condition at r = 0, the PMS length scale has to be smaller than r. 

Table [I] shows the errors e„ = E^ pprox — jrjexact f Qr ^ e Sj_p, and d states of the Morse potential with n = and 
n = 5. It compares present results with those of Baye et al y|. The last column of Table [I] displays the PMS optimal 
values of the grid spacing. 

Fig. [S] shows the local error t]n{x) = \ip( N \x) — ?/>( 80 ) (x)\ for the ground state of the Morse potential and approx- 
imations TV = 20 and N = 40, where we have chosen r = 3. We assume that the approximation of order N = 80 is 
sufficiently close to the exact wavefunction. Our results are more accurate than those of Baye et al Q who essentially 
considered the average of rj{x) over a chosen region. The sharp drop of the error beyond a certain value of r, clearly 
noticeable in Fig. |H] is due to the fact that our LSF reproduce the wave function only in a finite region outside which 
r}{x) is merely the value of ijj^ S0 \x)- It is worth noticing that t]n(%) is almost uniform in the region covered by the 
LSF. 
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FIG. 8: Error t)n(x) = \ip^ '(x) — ip^ as a function of x for the ground state of the Morse potential (color online). 



We have also applied our method to the one-dimensional Morse potential considered by Wei [l5j : 



V{x) = D [e- 2ax - 2e- ax + l] , 



(49) 



where — oo < x < oo, D = 0.0224, a = 0.9374, m = 119406 and fi = 1. This problem can be solved exactly and the 
energies are given by 



E n = Hlo 



1 huj 



n H n H — 

2 4D V 2 



1 



(50) 



where u> = y2D/m a. As before we can improve our results by conveniently shifting the potential on the x-axis 
from V(x) to V(x + x). When N = 20 and assuming x = we have found that the first-excited-state eigenvalue 
is reproduced with an accuracy of about 8 x 10~ 14 , which is even better than the accuracy obtained by Wei with 
N = 64 3 . When N = 40 the error of our method is just 8.1 x 10 -21 for the ground state. 



V. CONCLUSIONS 



In this paper we have introduced a new class of sinc-like functions, which we name "little sine functions" (LSF), 
which share some properties of the usual sine functions, but arc defined on finite intervals. We have shown that the 
LSF collocation method provides accurate approximations for a wide class of problems when it is supplemented by 
the Principle of Minimal Sensitivity (PMS) that provides an optimal grid spacing. In particular, we have applied 
the LSF to the solution of the Schrodingcr equation with only bound states and with mixed discrete and continuous 
spectra. We have chosen benchmark models treated earlier by other authors and in all the cases we obtained more 
accurate results. It seems that the LSF collocation method is an interesting alternative algorithm for solving many 
mathematical and physical problems. 



3 Notice that Table I of Ref. Il ■"! omits the the ground state of the model. 
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